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1. Introduction 



Renormalization group functions play an important role in quantum field theory. They deter- 
mine the energy scale dependence of the parameters of the Lagrange density and thus are important 
tools to combine predictions of the theory from different energy regions. An important example in 
this respect is the inspection of the gauge coupling unification at high energies where precise exper- 
imental data at the electroweak scale combined with accurate calculations of the renormalization 
constants yields precise predictions. 

In this contribution we consider the renormalization group functions of the three gauge coup- 
lings in the Standard Model (SM) within the Modified Minimal Subtraction (MS) scheme. In this 
renormalization scheme the beta functions are independent of all mass scales present in the theory 
and it is thus relatively simple to solve the loop integrals. In fact, in the calculation presented here 
at most massless three-loop two-point functions have to be evaluated which are known for more 
than 30 years Actually the main difficulty of the calculation is the huge amount of contributing 
Feynman diagrams (about 10 6 ) which result from the large number of vertices and propagators. 
This requires an automated setup which not only generates and computes all Feynman diagrams 
but also provides the Feynman rules in an automated way. 

The results presented in these proceedings have been obtained in Refs. @ ||||. There have 
been a number of publications where the one- and two-loop expressions have been computed [Q, 
H|, q, 0, H[ 10, [ll], 12, |i3| ]. Also several three-loop results have been computed since the end of 



the seventies [14, 15, 16, 17, |T^, 19]. Four-loop corrections to beta functions are only known for 



QCD [Efl, 21] 



Let us in a first step define the couplings we want to consider. It is convenient to use instead of 
the fine structure constant 0!qed an d the weak mixing angle dw the gauge couplings in a SU(5)-like 
normalization given by 

3 cos 2 dw snr dw 

Note that it is straightforward to obtain the beta functions for 0!qed an d sin 2 dw once the j3« are 
know (see, e.g., Ref. [^]). 



2. Calculation 



We define the beta functions via 

M 2 ^^=A({«>)},e), (2.1) 

where i = 1,2,3 labels the three gauge couplings. The index j runs over all couplings in the SM, 
i.e., the gauge, Yukawa and Higgs boson self couplings. We furthermore have e = (4 — d)/2 where 
d is the space-time dimension used for the evaluation of the momentum integrals. 

The functions /},•({ G5/(/i)},£) are conveniently computed from the renormalization constants 
relating the bare and renormalized gauge couplings via 

a, bare = At 2e Z„, ({«,}, £)«,-. (2.2) 
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Inserting this equation into (lO) and exploiting the fact that af aie does not depend on jx leads to 



£ Ck Oj_ y dZ Uj 
% Z a . f* . d a 



1 + 

Z a . dcij 



which constitutes the master formula for our enterprise. From this equation it is clear that it is 
sufficient to compute the renormalization constants Z Uj in the MS scheme in order to obtain jS,. In 
fact, we have to compute Z U[ , Z ai and Z U3 to three-loop order and the renormalization constants 
for the Yukawa coupling to one-loop order. As far as the Higgs boson self coupling is concerned 
it is sufficient to have the leading term proportional to e of the corresponding beta function. The 
discussion in the following is centered around the three-loop calculation of the gauge coupling 
renormalization constants. 

The procedure for the calculation of Z Uj (i = 1,2,3) is as follows: (i) choose a vertex which 
contains the coupling a,-; (ii) compute the renormalization constant of that vertex, Z vrtx ; (iii) com- 
pute the wave function renormalization constant of the external particles, Zt iW f ; and (iv) combine 
them according to Z«, = (Z vrtx ) 2 / '(Ht z *,wf) ■ 

We have used two independent approaches to compute the beta functions. The first one is 
based on the formulation of the SM using Lorenz gauge in the unbroken phase, i.e., all particles are 
still massless. Since the beta functions are mass independent this setup is quite convenient as the 
structure is simpler than after spontaneous symmetry breaking. The gauge bosons in this approach 
are the B and W bosons and the gluons. 

For the computation of Z Ul and Z aj it is convenient to choose the gauge boson ghost vertex 
since in that case the number of contributing diagrams is smallest. In order to have a cross check 
we have chosen the triple gauge boson vertex as a second option. For Z« 2 also the + 0~W3 vertex 
as a third alternative has been used. Note that due to the Ward identity Z ai is solely obtained from 
the B boson two-point function. Sample Feynman diagrams for all two- and three-point Green's 
functions can be found in Fig. [jj. 

In the second approach we have used the SM Lagrangian in the background field gauge (BFG) 



as a starting point (see, e.g., Ref. Q22|]). The advantage of this method is that only the gauge boson 
two-point functions and thus less different Green's functions have to be considered. On the other 
hand more vertices are present and thus more Feynman diagrams contribute to the background 
gauge boson propagator as compared to the corresponding quantum version. We use the BFG in 
the broken phase of the SM. However, here we can also set the masses of all particles to zero. 
Sample Feynman diagrams for the gauge boson two-point functions are shown in the first and 
second row of Fig. [I] where the external lines correspond to background gluon, photon and W or Z 
bosons. 

An important issue for all loop calculations involving electroweak gauge bosons is the treat- 
ment of y 5 in d ^ 4 dimensions. In our calculation we on purpose do not use Green's functions 
involving external fermions. For all other two- and three-point functions we could show that a naive 
treatment of y 5 leads to the correct result pl|. Note that Green's functions with external fermions 
are unavoidable for the calculation of the beta functions for the Yukawa couplings. At two-loop 
order it is again possible to use a naive version of 75. Beyond two loops, however, a more careful 



treatment is required (see, e.g., Ref. [23]). 
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Figure 1: Sample Feynman diagrams contributing to the Green's functions which have been used for the 
calculation of the renormalization constants of the gauge couplings. Solid, dashed, dotted, curly and wavy 
lines correspond to fermions, Higgs bosons, ghosts, gluons and electroweak gauge bosons, respectively. 
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We have performed several checks which convinced us from the correctness of our results. In 
brief they are given by: 



• comparison of the one- and two-loop results with the literature 



• comparison of partial three-loop results with the literature 



• computation of the beta function for the Higgs boson self-coupling to one-loop order and 
comparison with the literature 



• computation of the Yukawa beta functions for the top and bottom quark, and the tau lepton 
to two-loop order and comparison with the literature 



• computation of the BBB vertex to three-loop order; we checked that the sum of all 358 716 
diagrams gives zero 



• computation for general gauge parameters; we checked that they drop out in the final result 
for f3 functions 



• check that no infra-red divergences are present in the loop integrals 



As a last comment on the technical details of our calculation let us mention that we were able to 
obtain the final result for the gauge coupling beta functions for a general Yukawa structure involving 
all nine Yukawa couplings of the SM and the CKM matrix in the quark sector. Furthermore, it is 
straightforward to extend our final result for a SM with fourth family as will be explained in the 
next section. 



3. Results 



In this section we briefly discuss the final results for the gauge coupling beta functions. In 
order to show the structure of the analytical expressions we present the result for j8i which is given 
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(3.1) 



where the limit £ — > has been taken. For each gauge coupling beta function we have that the 
one-loop contribution of j8, is proportional to a 2 . Mixed contributions of order a 2 a 7 and afcCkCCi 
only appear at two and three loops, respectively, where a ; - are gauge or Yukawa and a* and a/ 
are gauge, Yukawa or Higgs boson self couplings. Note that the latter appears for the first time at 
three-loop order. 

The symbol no in Eq. (3J_) stands for the number of generations; in the SM we have hq = 3. 
The quantities trfi, tr T and trL incorporate the Yukawa couplings where we refer to Ref. [§] for 
details. In this contribution we only want to mention that the replacements 



trL" 



a" 



trf" 



a 



trS"->a£, trfs^a f a fo , 



(3.2) 



leads to the result where only the Yukawa coupling of the third generation is kept non-zero. 

In a similar way we can accommodate a fourth generation of fermions. In this case, the Yukawa 
matrices become 4x4 dimensional. If we assume that the fourth is much heavier and if we neglect 
all SM Yukawa interactions it contains a 3 x 3 zero matrix and we have 




with F = T,B,L, 



(3.3) 



where T and B stand for the up- and down-type heavy quarks, and L for the heavy charged and 

denotes the corresponding Yukawa couplings. Note that the 



aur. 



2 sin 2 ft 



neutral leptons, while ap 

contribution of a heavy neutrino is not contained in our formulae. 

Let us finally briefly discuss the numerical impact of the new three-loop corrections. In Fig. || 
we show the running of ai and a 2 from /I = Mz to the energy scales where these two coupling 
become equal. The dotted and dashed lines correspond to one- and two-loop running, respectively. 
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Figure 2: The running of the electroweak gauge couplings in the SM. The lines with positive slope corre- 
spond to 0Ci, the lines with negative slope to a-i. The dotted, dashed and solid lines correspond to one-, two- 
and three-loop precision, respectively. The bands around the three-loop curves visualize the experimental 
uncertainty. 

One observes a significant change of the curves, which is in particular much bigger than the experi- 
mental uncertainty indicated by the dashed band. Thus in case only one- and two-loop perturbative 
corrections are included the theory uncertainty is much bigger than the experimental one. This 
changes with the inclusion of the three-loop terms. The results are shown as solid lines which are 
close to the corresponding dashed curves. The effect is small, however, still of the order of the 
experimental uncertainty, in particular for 0C2- 

To conclude, the complete three-loop corrections to the gauge coupling beta functions have 
been computed [Q, |3]] using different methods and applying several cross checks. They constitute 
fundamental quantities of the SM and are important for high-precision experimental checks. 
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